Hybrid Quantum Gowdy Cosmology: Combining Loop and Fock Quantizations 



oo : 
o . 
o ■ 

(N ; 

3 : 

<: 
t> ■ 



cr 

i , 



(N 
> 
00 
On 
O 

O 

oo 
O 



X 



M. Martin-Benito 1 , L.J. Garay 1 ' 2 , and G.A. Mena Marugan 1 
1 Institute) de Estructura de la Materia, CSIC, Serrano 121, 28006 Madrid, Spain 
2 Departamento de Fisica Teorica II, Universidad Complutense de Madrid, 28040 Madrid, Spain 

We quantize an inhomogeneous cosmological model using techniques that include polymeric quan- 
tization. More explicitly, we construct well defined operators to represent the constraints and find 
the physical Hilbert space formed by their solutions, which reproduces the conventional Fock quan- 
tization for the inhomogeneities. The initial singularity is resolved in this inhomogeneous model in 
an extremely simple way and without imposing special boundary conditions, thus ensuring the ro- 
bustness and generality of this resolution. Furthermore this quantization constitutes a well founded 
step towards the extraction of physical results and consequences from loop quantum cosmology, 
given the central role of the inhomogeneities in modern cosmology. 

PACS numbers: 04.60.Pp, 98.80.Qc, 04.62. +v 



Loop quantum cosmology (LQC) |l[ has recently un- 
dergone important developments. The study of the phys- 
ical consequences of loop quantum gravity (LQG) [2| for 
cosmology has motivated the analysis of symmetry re- 
duced models that can be fully quantized by techniques 
that closely resemble those of LQG Q . Various homoge- 
neous and isotropic models have been successfully quan- 
tized (3). Other less symmetric systems have been also 
analyzed 0,0,0] with partially satisfactory results. 

Any realistic attempt to extract physical results from 
LQC must consider inhomogeneities, which play a cen- 
tral role in current cosmology. As a most suitable arena 
to progress in this direction, we will study the quantiza- 
tion of the linearly polarized Gowdy T 3 model [8] . This 
family of spacetimes is highly symmetric, but preserves 
one important feature of the full theory: it has an infinite 
number of degrees of freedom. Its classical solutions are 
indeed well known and generically present an initial sin- 
gularity [9j. Besides, its quantization has been addressed 
since long time ago by means of standard non-polymeric 
techniques [l(| (see also [ll| for a related system, namely 
the cylindrical waves). More recently, a rigorous Fock 
quantization of this model has been accomplished, which 
is shown to be essentially unique Nonetheless, even 
this Fock quantization (performed after deparametrizing 
the system) fails to resolve the cosmological singularity. 

To overcome this problem, we will carry out a hybrid 
quantization, that combines the polymeric quantization 
characteristic of LQC applied to the homogeneous solu- 
tions (which describe Bianchi I spacetimes) with Fock 
quantization of the inhomogeneities. This approach in- 
vestigates the effects on quantum geometry underlying 
LQC only on the homogeneous sector, while disregards 
the discreteness of the geometry encoded by the inho- 
mogeneities. A most natural treatment for the inho- 
mogeneities is then the Fock quantization. Indeed, one 
would expect that a quantum field theory for the inho- 
mogeneities, which can be regarded as a field living on a 
homogeneous (Bianchi I) background, be approximately 
valid on the polymerically quantized background. An 



additional motivation for this approach is the possibility 
that the inhomogeneities can be finally considered as per- 
turbations. We will represent the quantum constraints 
of this model, find their solutions, construct the physi- 
cal Hilbert space, and recover the conventional quantum 
field theory. To our knowledge this is the first time that 
this task has been carried out to this level of completion 
in a system with an infinite number of physical degrees 
of freedom. 

With this aim, we will reduce the classical phase space 
imposing gauge fixing conditions but, unlike in [12], we 
will not deparametrize the system completely, so that 
we will be left with a remanent Hamiltonian constraint, 
namely the spatial average of the full Hamiltonian con- 
straint. Its quantum counterpart is not entirely derived 
from LQG, since the inhomogeneities are not quantized 
a la loop, but it retains interesting polymeric features. 
In particular, we will face the problem of polymerically 
quantizing the internal time. As we will see this may 
affect the conventional concept of evolution. 

We will also analyze the role played by the inhomo- 
geneities in the quantum theory and thus the robustness 
of the results obtained in homogeneous LQC. In fact, as 
we have commented, a fundamental motivation for this 
hybrid quantization is the resolution of the cosmological 
singularity. The first studies in homogeneous LQC al- 
ready indicate that the section of constant internal time 
corresponding to the initial singularity can be avoided in 
the quantum theory 0, [l3| . In this letter we will show 
that the loop quantization of the homogeneous degrees 
of freedom, which correspond to global variables, suffices 
to "cure" the singularity, without the need to resort to 
the BKL (Belinsky, Khalatnikov, and Lifshitz) conjecture 
[13 ], Therefore we do not have to consider all the points 
individually. We will see that the initial singularity dis- 
appears from our quantum theory. Actually, nontrivial 
physical states contain no contribution corresponding to 
the singular initial section. Furthermore, they do not 
even evolve to another branch of the universe with oppo- 
site (triad) orientation, and can thus be viewed as arising 
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from "nothing" , without imposing any particular bound- 
ary condition. 

The considered Gowdy model is a globally hyperbolic 
vacuum spacetime, with spatial sections homeomorphic 
to a 3-torus and two hypersurface-orthogonal spatial 
Killing fields [§]. It can be described in terms of coor- 
dinates {t, 9, a, 5} adapted to the symmetries, where d a 
and d$ are the two Killing fields, with a, 5 6 S . The in- 
duced spatial 3-metric, the densitized lapse function N, 
and the shift vector N l , with i 6 {8, a, 8}, depend on t 
and 6 G S , and are then expandable in Fourier series. 

Suitable gauge fixing conditions for the diffeomor- 
phisms generated by the two Killing fields, together with 
the linear polarization condition, lead to a diagonal 3- 
metric. It can be characterized by three fields that de- 
scribe the norm of one of the Killing vectors, the area 
of the isometry group orbits, and the scale factor of the 
metric induced on the set of group orbits. We further 
impose that the generator of the conformal transforma- 
tions of this latter metric and the area of the isometry 
group orbits be homogeneous functions. These condi- 
tions fix the gauge freedom associated to the nonzero 
Fourier modes of the 0-momentum constraint and of the 
densitized Hamiltonian constraint, and imply that the 
functions N 6 and N are homogeneous [Hj]. The classi- 
cally reduced system is then described by three canonical 
pairs of variables, which correspond to the homogeneous 
degrees of freedom of a Bianchi I spacetime with T 3 - 
topology, and by the infinite number of nonzero (inho- 
mogeneous) modes of the field unaffected by the gauge 
fixing, together with their corresponding conjugate mo- 
menta (we adopt the same parametrization as in [lij]). 
We will call these the homogeneous and inhomogeneous 
sectors respectively. There remain two constraints: the 
zero mode of the ^-momentum constraint and the zero 
mode of the densitized Hamiltonian constraint, which 
generate S 1 and time translations. 

The homogeneous sector of this reduced Gowdy model 
is described in Ashtekar variables. We adopt a diagonal 
gauge for these variables, corresponding to Bianchi I de- 
grees of freedom. Since they can be defined in such a way 
that they do not depend on the selected fiducial metric 
Q , we choose the Euclidean one for simplicity. The nat- 
ural coordinate cell of the system is the T 3 -cell. Then the 
nontrivial components of the SU(2) gravitational connec- 
tion and of the densitized triad are c 1 /(2tt) and Pi/(47r 2 ) 
respectively, with {c l ,pj} = 8irGj5j. Here, 7 is the Im- 
mirzi parameter. We note that pg is proportional to the 
time variable chosen in (l2T | to deparametrize the system. 

The kincmatical Hilbert space for the homogeneous 
sector is constructed adapting LQG procedures to sym- 
metry reduced systems 0, [BJ . The elementary homoge- 
neous variables are holonomies along edges of oriented 
coordinate length 2-K\ii in the i direction, ^ being a real 
number, and triad fluxes through rectangles orthogonal 
to those directions. The configuration algebra, denoted 



by Cylg, is the algebra of almost periodic functions gen- 
erated by the exponentials Af Mj (o 3 ) = exp(ifj,jd' /2) [15[ , 
which in the Dirac ket notation become \/J,j). The ho- 
mogeneous sector of the kinematical Hilbert space is the 
Cauchy completion of Cyl s with respect to the discrete 
inner product (MilMi) = <^W'' f° r eacri direction. The ba- 
sis states are eigenstates of the momentum operator 
Pi, while the action of .A/O shifts their label \Xi. 

Remember that in LQG physical areas are discrete and 
have a minimum nonzero eigenvalue A. This feature has 
an imprint in LQC. The coordinate length of the holon- 
omy along each edg e exhibits a state-dependent mini- 
mum value 2np,i [16[. Here, we will adopt the (operator) 
relation p>f\pi\ = A, as proposed in Q (for additional 
details on this relation, see [17|). As a result, the opera- 
tor generates a state-dependent minimum shift. It is 
convenient to rename the states by reparametrizing the 
labels fii so that this minimum shift becomes uniform. 
Carrying out this affine reparametrization 0, [B[ one in- 
troduces, for each direction, a new label Vi(/ii) such that 
tffiifri) = l«< + 1} and pi\vi) = 3 1 / 3 Asgn(w i )|w i | 2 / 3 |u i ). 
The operators acting on the homogeneous sector can then 
be written in terms of these elementary ones, which are 
densely defined in the domain spanned by the states \vi). 

For the inhomogeneous modes of the only field present 
in the theory, we adopt the Fock quantization of [12j . 
We represent this field in terms of the basis of n-particle 
states naturally associated to a free massless scalar field 
03, \{n m }) ■= |...,n_ m , ...,n m , ...), where m ^ and 
n m < 00 is the occupation number of the m-th mode. In 
these states, only a finite set of these occupation numbers 
differ from zero. The proper Fock subspace annihilated 
by the constraint that generates S^-translations will be 
called J- p . A basis for it is provided by the n-particle 
states that verify the condition Ylm>o m ( n m — n-m) = 0. 

In order to quantize the densitized Hamiltonian con- 
straint of the model, let us first focus on Bianchi I. 

LQG techniques provide the (nondensitized) quantum 
Hamiltonian constraint Cbi for Bianchi I in terms of the 
elementary operators pi and [1, [13, EH • We adopt 
a factor ordering such that this constraint is a sum of 
terms factorized in each direction, each of these factors 
being symmetric. Besides, our prescriptions ensure that 
Cbi annihilates the proper subspace on which any of the 

Taking into account that 
we will call it the 



operators pi vanishes [17 
the volume operator is V 
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subspace of zero- volume states. Remarkably, this sub- 
space decouples when the constraint is imposed because, 
in addition, its complement is invariant under the action 
of Cbi- Therefore, we will restrict to this complement in 
the following. 

We are now ready to write down the densitized quan- 
tum Hamiltonian constraint Cbi for Bianchi I. LQG pro- 
cedures give the regularized inverse volume operator 1 /V. 
It only annihilates the zero-volume states, which are re- 
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moved from our theory. Consequently, its inverse [1/V] 
is well denned. In principle, Bianchi I physical states (i.e. 
those annihilated by Cbi) will not be normalizable in the 
kinematical Hilbert space. These states, denoted by (tp\, 
live in a larger space on which we can define 

Cbi = [l/V}- 1/2 C B i[l/Vr 1/2 - 

This operator annihilates the transformed physical states 
(-01, related to (-01 by the bijection $\ = (ifj\ [ljv] 1/2 . 

In the full Gowdy model, once the homogeneity of the 
densitized lapse function is taken into account, the den- 
sitization of the Hamiltonian constraint is carried out by 
means of the transformation explained above. Again, 
states on the kernel of any of the p^s are decoupled. The 
densitized Hamiltonian constraint Cq is the sum of two 
terms: the already quantized Bianchi I part and the part 
that involves the inhomogeneities. The latter depends on 
the absolute value of our internal time variable, \pg\- In 
particular, there is a term quadratic in the inverse of its 
square root, for which LQC provides the corresponding 
regularized operator (which acts diagonally on the states 
\vg)). In addition, this part of the constraint also depends 
on the connections, which are not well defined operators 
in LQC. This difficulty can be overcome if we notice that 
the only dependence on the connections occurs through 
factors of the form c l pi. Indeed, if we compare Cbi with 
its classical counterpart, it is possible to identify a natural 
quantum prescription to represent dpi . The correspond- 
ing symmetric operator, 0j, can be written in terms of 
the operators \pi\, sgn(pj), and jV±2pj- It can be shown 
to be essentially self-adjoint. Its action on the states \vi) 
is of the form 0, EH 

6iK) = -i-K^l [f+(vi)\vi + 2) - f-(vi)\vi - 2)] , 

where l p = \J Gh is the Planck length and the functions 
f+(vi) and f-{vi) possess the relevant property that they 
vanish in the intervals [—2,0] and [0,2], respectively. 
The explicit expression of the constraint Cg is [HI 
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is the Hamiltonian of a mass- 



Here, fl| = £„. 
less free scalar field, and a m and aj^ are the annihilation 
and creation operators of particles in the m-th mode cor- 
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responding to that field, such that [d„ 
the above formula for Cg, the first line is the densitized 
Hamiltonian constraint Cbi of Bianchi I. 

Since 0i is a difference operator and actually it does 
not relate states IvA with Vi > to states with v; < 0, 



the label Vi can be restricted to any of the semilattices 
C% = {±(e l + 2k), k e N}, where S (0,2]. Semilat- 
tices corresponding to different E{ or different signs do 
not get mixed under the action of the constraint. In this 
sense the homogeneous sector is superselected. We can 
then focus our study on the kinematical Hilbert space 
(®iH ei ) <S> !F P , where H ei is the completion of the span 
of states \vi) with Vi in the specific semilattice £+. 

Note also that CT and 05 are Dirac observables. A 
careful analysis, complemented with numerical evidence, 
supports that the spectrum of 0^ is absolutely continu- 
ous and coincides with the real line [13] ■ Therefore its 
generalized eigenstates |Ai) are delta-normalizable with 
respect to the Lebesgue measure d\. 

The constraint Cg is a well defined operator in the do- 
main spanned by {\vg)®\Vcr)®\v5)®\{n m })\Vi ^ 0} [20] ]. 
This result shows the viability of combining the poly- 
meric and Fock quantizations. Furthermore, an analysis 
of the associated deficiency index equation indicates that 
it possesses no normalizable solutions and therefore that 
the constraint operator is essentially self-adjoint. Thus, 
group averaging techniques can in principle be available 
to determine the space of solutions to the Hamiltonian 
constraint. 

Nonetheless, one can formally obtain the generic ex- 
pression of the solutions without the need to resort to 
group averaging, whose application is not immediate in 
this case. For convenience, we will use the (generalized) 
eigenstates of pg, Q a , and 0,5, instead of those of the p^'s, 
to decompose the homogeneous sector. Then, if we for- 
mally expand the states (ip\, we can represent them by 
the coefficients (ip\ a ,\ s (vg)\{n m }) . <dg is the only oper- 
ator in Cg that does not act diagonally on the homoge- 
neous sector. Consequently, the Hamiltonian constraint 
becomes a recurrence equation that relates the coeffi- 
cients in the three consecutive sections (vg — 2), vg, and 
(vg + 2). However the coefficient evaluated at (vg — 2) ap- 
pears multiplied by the factor f_(vg). Since f~{£g) = 0, 
the recurrence equation turns into a consistency relation 
for the coefficients in the two first sections eg and (eg + 2). 
As a result any coefficient (ip\ a } \ s (vg)\{n m }) can be ex- 
pressed in terms of the initial ones (i(jx„,\ s (eg)\{rifh})- 
The formal solution of the constraint is, obviating the 
A-dependence and the projection on n-particle states, 



(4>(es +2k)\ = $(eg)\ [ II F ( £t > + 2r * + 2 ) 

0(k) {r,} 

xv\ i Y[H i (e e + 2s 



where 0(k) is the set of possible ways to move from 
to k by jumps of one or two steps and, for each element 
of O(k), {sj} is the set of integers followed by a jump of 
one step, while {r^} are those followed by a jump of two 
steps [20j. Besides, F(vg) = f-(vg)/ f+(vg), the symbol 
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V denotes path ordering, and H^(vg) has the form: 

H s (ve) = yiM + yi(ve)H^ + y 3 (v e )Hl, 

Uk(,vg) being A-dependent functions whose detailed ex- 
pression is irrelevant for our discussion. 

As we have seen, {tp\ at \ s {sg)\ completely determines 
the solution, so that we can identify the latter with the 
corresponding data on the initial section vg = eg. To 
endow the space of these solutions with a Hilbert struc- 
ture, we choose a complete set of real classical observ- 
ables and require that their quantum counterparts be 
self-adjoint operators. At this stage it is easier to relax 
the S' 1 -gauge symmetry, so that we can choose the op- 
erators that represent the Fourier sine and cosine coeffi- 
cients of the nonzero modes as the basic observables for 
the inhomogeneous sector. For the homogeneous sector, 
we can select as a complete set of observables (in the A- 
representation) the operators of multiplication by A CT and 
Xs (i.e., CT and 0a), together with the respective deriva- 
tives evaluated on the considered section, —id\ a \ Eg and 
— id\ s \e e - I n order to obtain the physical Hilbert space 
we first demand reality conditions on all these observ- 
ables and then restore the ^-symmetry. The final result 
is that, for generic initial data (ip\ lrt \ s (sg)\, the physical 
Hilbert space that one obtains in this way is unitarily 
equivalent to L 2 (R 2 , dXvdXs) <£> T v . In particular, let us 
emp hasize that we then recover the Fock quantization of 
[12| for the physical inhomogeneous sector. 

Several comments are in order once we have completed 
the quantization of this Gowdy model. 

We remark that the general expression of the solutions 
is only formal. A more rigorous interpretation in terms 
of structures associated to the kinematical Hilbert space 
is obstructed by the production of infinite particle pairs 
caused by the operator H^ nt and by the t^-dependence of 
the inhomogeneous part of the constraint, which prevents 
its diagonalization in terms of particle states with a fixed 
mass for all the fg-sections. Nevertheless we can in fact 
make full sense of the solutions in certain situations. 

This is the case if one introduces a cut-off for the 
wavenumber of the field modes, what can be understood 
as a further symmetry reduction of the Gowdy model. 
Solutions to the constraint are then in the algebraic dual 
of the dense set of n-particle states on each vg-section. 
Actually, the mentioned physical Hilbert space for the full 
Gowdy model can then be attained in the limit in which 
the cut-off is removed. A related possibility that we will 
explore in a future research is to introduce a dynamical 
cut-off, dependent on the value of vg, as an effective way 
to account for the minimum coordinate length fig. 

Other situations in which the structure of the solutions 
is well under control is when one regards the inhomoge- 
neous part of the Hamiltonian constraint as a perturba- 
tion around Bianchi I and truncates the solutions at sec- 
ond perturbative order in it or, alternatively, if one views 



as a perturbation only the term that contains iT?L and 
makes the truncation at first order. Then the solutions 
can be defined again, up to the analyzed perturbative 
order, in the dual of the n-particle states for each in- 
sertion and the physical Hilbert space that one obtains 
is just L 2 (M 2 , dX a dX s ) T p . 

Note that in the above cases we can adopt the picture 
that solutions are initial data on the section vg = eg which 
evolve to new data on any other section vg — eg + 2k. 
Identifying sections, we arrive in principle at a concept of 
evolution for physical states which would not be unitary 
(i.e. the norm would not be preserved) as a consequence 
of the polymeric quantization of the internal time. 

Furthermore, the loop quantization carried out on the 
homogeneous sector allows us to get insight about the 
resolution of the cosmological singularity. The classical 
Gowdy model generically presents an initial singularity 
where some of the variables pi vanish. Since zero-volume 
states are removed in the quantum theory, the kernel of 
any of the operators pi is empty. Therefore the singular- 
ity disappears already at the kinematical level, and, at 
least in this sense, is resolved. Note that, in a standard 
quantization of the homogeneous sector, the operators pi 
have instead a continuous spectrum. Thus, their kernel 
is not a proper subspace which could be decoupled in or- 
der to remove the singularity in the kinematical Hilbert 
space. Moreover, the singularity persists then also in the 
physical Hilbert space. This fact can be easily realized 
in the completely deparametrized system quantized a la 
Fock [l2j , since the domain of definition of the time pa- 
rameter plays there the role of the spectrum of \p$ | quan- 
tized in the standard way. Curvature invariants become 
explicitly time-dependent quantum observables that in- 
deed diverge in the limit when the time parameter van- 
ishes [§] . 

To conclude, we have successfully quantized this inho- 
mogeneous Gowdy model by combining loop and Fock 
quantizations. The attained results prove that our hy- 
brid quantization is feasible and well posed. More pre- 
cisely, we have represented rigorously the quantum con- 
straints and explicitly found their solutions, which have 
been endowed with a Hilbert structure. The physical 
Hilbert space thus obtained reproduces the conventional 
Fock quantization of the inhomogeneous modes. As far as 
we know, this is the first quantization carried out to this 
degree of completion in an inhomogeneous cosmological 
system using polymeric techniques. The loop quantiza- 
tion of the homogeneous sector proves enough to resolve 
the classical singularity, inasmuch as it disappears from 
the quantum theory. In fact, it is not only avoided in 
physical states, but they never cross it. The quantum 
states of our Gowdy universe start in a single nonsingular 
section, without requiring any boundary conditions. This 
result reinforces the robustness of the big-bang avoidance 
in LQC and paves the road for the extension of the sin- 
gularity resolution to inhomogeneous cases. 



5 



The authors are very grateful to J.M. Velhinho and 
T. Pawlowski. This work was supported by the Spanish 
Grants FIS2005-05736-C03-02 (its continuation FIS2008- 
06078-C03-03), FIS2006-26387-E, and CSD2007-00042 
(CPAN); and M.M-B. by CSIC and the European Social 
Fund under the grant I3P-BPD2006. 



[1] M. Bojowald, Living Rev. Rel. 11, 4 (2008). 

[2] T. Thiemann, Modern Canonical Quantum General Rel- 
ativity (CUP, Cambridge, UK, 2007). 

[3] A. Ashtekar, M. Bojowald, and J. Lewandowski, Adv. 
Theor. Math. Phys. 7, 233 (2003). 

[4] See, e.g., A. Ashtekar, T. Pawlowski, and P. Singh, Phys. 
Rev. Lett. 96, 141301 (2006); Phys. Rev. D 73, 124038 
(2006); 74, 84003 (2006). 

[5] D.W. Chiou, Phys. Rev. D 75, 24029 (2007). 

[6] D.W. Chiou, Phys. Rev. D 76, 124037 (2007). 

[7] K. Banerjee and G. Date, Class. Quantum Grav. 25, 
145004 (2008). 

[8] R.H. Gowdy, Ann. Phys. 83, 203 (1974). 

[9] V. Moncrief, Phys. Rev. D 23, 312 (1981). 
[10] See, e.g., C.W. Misner, Phys. Rev. D 8, 3271 (1973); 
B.K. Berger, Ann. Phys. 83, 458 (1974); Phys. Rev. D 
11, 2770 (1975); Ann. Phys. 156, 155 (1984); GA. Mena 
Marugan, Phys. Rev. D 56, 908 (1997); M. Pierri, Int. J. 



Mod. Phys. Phys. D 11, 135 (2002). 
[11] See, e.g., K. Kuchar, Phys. Rev. D 4, 955 (1971); A. 

Ashtekar and M. Pierri, J. Math. Phys. 37, 6250 (1996); 

J.F. Barbero G, GA. Marugan, and E.J.S. Villasenor, 

Phys. Rev. D 67, 124006 (2003). 
[12] A. Corichi, J. Cortez, and GA. Mena Marugan, Phys. 

Rev. D 73, 41502 (2006); 73, 84020 (2006); A. Corichi, J. 

Cortez, GA. Mena Marugan, and J.M. Velhinho, Class. 

Quantum Grav. 23, 6301 (2006); J. Cortez, GA. Mena 

Marugan, and J.M. Velhinho, Phys. Rev. D 75, 84027 

(2007). 

[13] M. Bojowald, Phys. Rev. Lett. 86, 5227 (2001). 

[14] V.A. Belinsky, I.M. Khalatnikov, and E.M. Lifshitz, Adv. 
Phys. 31, 639 (1982). 

[15] We will not use the Einstein summation convention. 

[16] There exist different proposals to define this minimum 
value, which are still under discussion. This issue has 
been considered in [|| and by A. Ashtekar (unpublished). 

[17] M. Martin-Benito, GA. Mena Marugan, and T. Paw- 
lowski, |arXEi0lQli3157| [Phys. Rev. D (to be published)]. 

[18] M. Martin-Benito, Cuantizacion del Modelo de Gowdy 
T 3 con Polarization Lineal (Master thesis, Univ. Com- 
plutense de Madrid, 2007). 

[19] Namely, the algebraic dual of the subset of Cyl s of 
nonzero-volume states. 

[20] L.J. Garay, M. Martin-Benito, GA. Mena Marugan, and 
J.M. Velhinho (unpublished). 



